Existence of Density Functionals for Excited States and Resonances 
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We show how every bound state of a finite system of identical fermions, whether a ground state 
or an excited one, defines a density functional. Degeneracies created by a symmetry group can 
be trivially lifted by a pseudo-Zeeman effect. When complex scaling can be used to regularize a 
resonance into a square integrable state, a DF also exists. 
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The aim of density functional (DF) theory is to con- 
tract a functional that provides the energy expectation 
value for a correlated many-body state as a function 
of the one-body density, such that minimization of the 
DF leads to the exact ground state (GS) density. Since 
the existence theorem proven for GSs by Hohenberg and 
Kohn (HK) Q, its extension by Mermin [2| to equilib- 
rium at finite temperatures, and the further development 
by Kohn and Sham (KS) [3| of an equivalent, effective, 
independent particle problem, a considerable amount of 
work has been dedicated to generalizations such as spin 
DFs |4(, functionals taking into account the symmetries 
of the Hamiltonian calculations of excited state den- 
sities la, |3 , treatments of degeneracies or symmetries of 
excited states 0, Q and quasiparticles ■ For the reader 
interested in an even more complete reading about both 
basic questions and applications, we refer to |llj-|lflj. 

DFs for resonant states have received much less atten- 
tion. We want to study this problem here. First we will 
address two related issues, namely that of a unified the- 
ory for ground and excited states and that of a theory 
for non degenerate and degenerate ones. A generalized 
existence theorem can be constructed by modifying the 
Hamiltonian in such a way that the spectrum is shuffled 
but the eigenstates are left unchanged, and by making 
a systematic use of the Legendre transform (LT) for a 
detailed analysis of the density. 

A reminder of the HK proof is useful here. Consider 
a finite number A of identical fermions, with at and ap 
their creation and annihilation operators at position r, 
and the physical Hamiltonian, H = T + V + U, where 

T = Y,i=iU, V = Et>j=i v ii and U = Yn=x u i are 
the kinetic, two-body interaction and one-body poten- 
tial energies, respectively. For simplicity, we consider 
such fermions as spinless and isospinless and work at 
zero temperature. Both v and u may be either local 
or non local. Next, embed the system into an addi- 
tional one-body, external field, W = J2t=i w ii to observe 
its (non linear!) response. The Hamiltonian becomes 



K = H + W. It is understood that w is local, (f\w\r') = 
w(r) S(f — r'), although a DF theory with non-local po- 
tentials exists j2jj. The usual Rayleigh-Ritz variational 
principle, where \ip) is just an A-particle, antisymmetric, 
square normalized, otherwise unrestricted wave function, 
applied to Fm = min^ F, with F = (?/>|K|t/;), generates 
ipmin, the exact GS of K, with the exact eigenvalue Fm- 
The minimum is assumed to be non degenerate, smooth, 
reached. Clearly, ipmin and Fm are parametrized by 
w. An infinitesimal variation 8w triggers an infinitesi- 
mal displacement Stpmin, with 5F M = (ipmin\SW\ip m i n ) . 
There is no first order contribution from 5ip m i n . De- 
fine the one-body density matrix in coordinate repre- 
sentation, n(f,r') = ("0mm | at a? / 1 V'rmn}- Its diago- 
nal, p{f) — n(f,r), is the usual density deduced from 
\^min\ 2 by integrating out all particles but one. Since 
SFm = J dr p{r) 5w{r), then SFm /Sw(r) = p(r). Freeze 
t, v and u and consider Fm as a functional of w alone. 
The HK process then consists in a Legendre transform 
of Fm, based upon this essential result, SFm/Sw = p. 
This LT involves two steps: i) subtract from Fm the 
functional product of w and SFm/Sw, i.e. the integral 
/ dfw{r) p(r), leaving T M = (ipmin |H|^ min ); then ii) set 
p, the "conjugate variable of w" , as the primary variable 
rather than w; hence see Tm as a functional of p. Step 
ii) is made possible by the one-to-one (1 <-> 1) map be- 
tween w and p, under precautions such as the exclusion 
of trivial variations 5w that modify w by a constant only, 
see for instance [lfi| and jlflj . The 1 <-> 1 map is proven 
by the usual argument ad absurdum : if distinct poten- 
tials w and w' generated VVmn and 4>min (distinct!) with 
the same p, then two contradictory, strict inequalities 
would occur, J df[w(r) — w'(f) ] p(f) < Fm — F M , and, 
/ dr[w(f) — w'(r)] p{r) > Fm — F M . An inverse LT re- 
turns from Tm to Fm, because STm/Sp — —w. Finally, 
the GS eigenvalue Eq of H obtains as Eq = min^ Tm [p] ', 
the GS wave function ipQ of H is the wave function i/Wn 
when w vanishes; that density providing the minimum of 
Tm is the density of ipo . 
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Consider now any excited bound eigenstate ip n of H, 
with its eigenvalue E n . Then, trivially, ip n is a GS of 
the semipositive definite operator (H — E n ) 2 . Since E n 
is not known a priori, consider rather an approximate 
value E n , obtained by any usual technique (configuration 
mixing, generator coordinates, etc.) and assume that E n 
is closer to E n than to any other eigenvalue E p . Then 
ip n is a GS of (H — E n ) 2 . The possible degeneracy de- 
gree of this GS is the same whether one considers H or 
(H - E n ) 2 . Introduce now K = (H E n f + W. If 
there is no degeneracy of either if> n or its continuation 
as a functional of w, then the HK argument holds as 
well for K as it does for K. Hence a trivial existence 
proof for a DF concerning ip n . But most often, ip n be- 
longs to a degenerate multiplet. Degeneracies are almost 
always due to an explicitly known symmetry group of 
H. Notice however that the external potential w does 
not need to show the same symmetry; hence, in general 
for K, there is no degeneracy of its GS; a unique ipmin 
emerges to minimize the expectation value of K. How- 
ever, for that subset of zero measure in the space of po- 
tentials where w shows the symmetry responsible for the 
degeneracy, and in particular for the limit w — > 0, pre- 
cautions are necessary. Consider therefore an (or several) 
additional label(s) g sorting out the members ip n g of the 
multiplet corresponding to that eigenvalue (E n — E n ) 2 
of (H — E n ) 2 . There is always an operator G related to 
the symmetry group, or a chain of operators Gj in the 
reduction of the group by a chain of subgroups, which 
commute with H and can be chosen to define g. For sim- 
plicity, assume that one needs to consider one G only. 
Then define g as an eigenvalue of G and assume, obvi- 
ously, that the spectrum of G is not degenerate, to avoid 
a reduction chain of subgroups. It is obvious that, given 
some positive constant C, and given any chosen 7 among 
the values of g, there is no degeneracy for the GS of 
(H — E n ) 2 + C (G — 7) 2 . Nor is there a degeneracy of the 
GS of K = (H-E n ) 2 + C (G- 7 ) 2 +W = K+C (G-7) 2 , 
even if w has the symmetry. When several labels be- 
come necessary with a subgroup chain reduction, it is 
trivial to use a sum • Cj (Gj — jj) 2 of "pusher" terms. 
Finally a DF results, now from the HK argument with 
K. We stress here that pusher terms, because they com- 
mute with H, do not change the eigenstates of cither 
H nor (H — E n ) 2 . Only their eigenvalues are sorted 
out and reorganized. Note that the pusher expecta- 
tion value vanishes for ip ni . Naturally, when w is fi- 
nite, eigenstates of K differ from those of K, but what 
counts is the information given by the DF when w van- 
ishes. A simplification, avoiding cumbersome square op- 
erators H 2 , is worth noticing. Consider the operator, 
K = H + C (G - 7) 2 + W. At the limit where w van- 
ishes, there is always a choice of a positive constant C 
which makes the lowest state with quantum number 7 
become the GS. This leads to a more restricted density 



functional that is of interest for the study of an yrast line. 

That DF, Tm[p\: based upon K, provides the ex- 
pectation value, T M \p\ = (VWi| [ (H - E n ) 2 + C (G - 
7) 2 ] \ipm,i n ), where ipmin, square normalized to unity, is 
also constrained by the facts that (VWn| tvoHVWn) = 
p(f) and K\ip m in) = e\ip mm ) for the eigenvalue e = 
Fm- It may be interesting to find a DF that pro- 
vides the expectation value of H itself. This can be 
done by taking the derivative of J-m[p] with respect 
to E n , at constant p. We suppose that this deriva- 
tive exists, which is the case for a discrete spectrum 
at least. With the notation \tp) — d\tp)/dE n , and us- 
ing the fact that (tp min | W|^ min ) + (Vw„|W|Vw) = 
fw(r) (dp(r)/dE n ^ df — 0, one can write: 

= 2 (Vwl (E n - H) |</w) + 

(lpmin\(s ~ W)|-0mm) + (lpmin\{£ ~ W)\lp m i n ) 

= 2 (VWi| [E n - |VWi) • (1) 
Therefore we can define a new DF, 



T D [p]=E n 



&Fm [p] 
2dE n 



(2) 



such that T D [p] = (tpmin \H-\ipmin) an d •? 7 D[Pn 7 ] = E n 
for the density p ni of the eigenstate i\) ni of H at energy 
E n . Furthermore one finds that ^^[pny] — 0, because 
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n = iJ J n~f ■ Hence the functional Td [p] is stationary 
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for ip. 

at the exact density p — p ni . It is not expected to be 
minimal at p ni , however, unless the resulting eigenstate 
corresponds to the absolute GS when w vanishes. 

Resonances may be defined as special eigenstates of 
H if one uses an argument a la Gamow, allowing some 
radial Jacobi coordinate r > to show a diverging, ex- 
ponential increase of the resonance wave function at in- 
finity of the form exp(ipr), where the channel momen- 
tum p is complex and 3p < 0. It is well known that 
those eigenvalues E n describing resonances are complex 
numbers, with 3-E„ < 0. There have been extensive dis- 
cussions in the literature about the physical, or lack of, 
meaning of such non normalizable wave functions and 
about the wave packets which might be used to replace 
them, HI El El 0. The point of view we adopt in 
this note is base d up on the Complex Scaling Method 
(CSM) jH Em ■ a modest modification of H trans- 

forms narrow resonances into square integrable states; 
then there is no difference between the diagonalization 
for bound states and that for resonances. The cost of the 
CSM, however, is a loss of hcrmiticity: the CSM Hamilto- 
nian H' is non hermitian, somewhat similar to an optical 
Hamiltonian jH HI |2| . 

Given the ket eigenstate equation, (H' — E n )\ip n ) = 0, 
where \ip n ) is now a square integrable resonance wave 
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function, we can consider the hermitian conjugate equa- 
tion, (^> n |(H't — E*) = 0. Clearly, tp n is a GS, as both a 
ket and a bra, of the hermitian and semipositive definite 
operator, Q exac t = (H't — E*) (H' — E n ), with eigenvalue 
0. Applying the same argument as before, but now to 
Qexact instead of (H — E n ) 2 , demonstrates the existence 
of a DF around the targeted resonant state. 

In practice we do not know E n exactly. Given a 
sufficiently close estimate E n of E n , an approximate 
GS eigenvalue \E n — E n \ 2 occurs for Q appr x = (H't — 
E*) (H' — E n ), at first order with respect to AQ = 
Qapprx — Qexact- Since ip n is not a ket eigenstate of 
H'" 1 * = H' — 2i SH', it is also perturbed at first order 
in AQ. Still one can copy the construction for Fd[p}, 
see Section 3, if one interprets the operator d/d-E* as 
d/ddi.E n + id/d^sE n . The resulting functional Toip] is 
linear in H'. For E n = E n the functional will be sta- 
tionary at the density of the exact resonant state. While 
provinding a proof of existence, the construction of the 
exact functional for H' requires the knowledge of the ex- 
act eigenvalue E n . This might be an inconvenient limita- 
tion but fortunately calculations of numbers such as E n 
are usually much easier and much more precise than cal- 
culations of wave functions ip n and/or of their densities. 

If the resonance has good quantum numbers (QN)s 
inducing degeneracies, the same pusher terms as those 
which have been discussed above can be added to create 
a unique GS, from the operator, Qexact + C (G — j) 2 . 
The HK argument, implemented with the full operator, 
K' = (H't - E*) (H' - E n ) + C (G - 7 ) 2 + W, then 
proves that DFs exist for those resonances regularized by 
the CSM. Notice, however, that a simplified theory, with 
an "yrast suited" operator K', linear with respect to H', 
is not available here, since the restoration of hermiticity 
forces a product H' 1 * H' upon our formalism. 

We now consider a special case of rather wide inter- 
est in nuclear and atomic physics, i) Good parity of 
eigenstates of H = T + V or H = H + U when 
u is restricted to be even, is assumed in the following. 
Hence our eigendensities, quadratic with respect to the 
states, have positive parities, ii) We also assume that 
the number of fcrmions is even, iii) The QNs in which 
we are interested in this Section are the integer angu- 
lar momentum L and magnetic label M of an eigenstate 
^ lm of H, where it is understood that the two-body v 
and one-body u interactions conserve angular momen- 
tum. When w is switched on and is not rotationally in- 
variant, eigenstates of K, K, or K may still tolerate such 
labels LM by continuity. First, consider w = 0. The 
density plm comes from the product ^ l * m ^lm, but it 
does not transform under rotations as an {LM} tensor. 
Rather, it is convenient to define "auxiliary densities", 
^Ao(r) = Em=-l(-) L - M (L - M L M | A 0) p LM (f), 
where (L — M L M | A 0) is a usual Clebsch-Gordan co- 
efficient. Each function <7 A o(r) now behaves under rota- 



tions as a {AO} tensor. It can therefore be written as the 
product of a spherical harmonic and a radial form factor, 
^Ao(r) = Y M) (f) r A (r) - y/(2X + l)/47r£ A (cos/3) r A (r), 
where C\ is a Lcgcndrc polynomial and the angle (3 is the 
usual polar angle, counted from the z-axis. Conversely, 

2L 

Plm(t) = ^(-) L - M (L-MLM|A0)r AO (r)r A (r). (3) 

A=0 

This provides a "Fourier analysis" of plm in angular 
space. The density is parametrized by scalar form fac- 
tors, t\. Since L is here an integer and furthermore pl-m 
and plm are equal, and since Clebsch-Gordan coefficients 
have the symmetry property (L M L' M' | A M") = 
(-) L+i '- A (L' M' L M | A M"), then necessarily t a = 
if A is odd. There are thus (L + 1) scalar functions, 
t , r 2 , t 2 l, to parametrize (L + 1) distinct densi- 
ties pLo, Pli, Pll- Because of the quadratic nature 
of the density observable, the even label A for angular 
"modulation" of p runs from zero to twice L, with a U 2L 
cut-off"; a signature, necessary if not sufficient, for an 
"L-density". Reinstate now w as the LT conjugate of 
Plm- It makes sense to study situations where w is re- 
stricted to expansions with (L + 1) arbitrary scalar form 
factors, w(r) = Y^lven A=o Y \o{r) w\(r). With inessential 
factors such as (-) L ~ M (L - M L M | A 0} omitted 
for simplicity in the following, every pair {rr\,rw\} is 
conjugate. An cigendensity of K, K, K may have an in- 
finite number of multipole form factors, but, with such 
restricted potentials w, only r ,r 2 , ...,t 2 l are chosen by 
the LT relating Tm and Fm- 

It can make even more sense to restrict w to one multi- 
pole only, w(r) = Y\o(f) w\(r), with A = 0, or 2, ... or 2L, 
to study each multipole of p separately. For simplicity we 
now use the easier version of the theory, with that oper- 
ator K which is suited to the yrast line. Add therefore 
to H a pusher term Zlm leaving intact the eigenstates, 
namely Z lm = B [t-t-L{L+l)] 2 +C (L 2 - M) 2 . Hence 
k LM x = H + Z LM + W A = T + V + U + Z LM + W A . 
Here the subscript A specifies that w is reduced to one 
multipole only. Then L is the total angular momentum 
operator and L z is its third component. This operator 
Z lm moves the eigenvalues of H so that the lowest eigen- 
state of H with quantum numbers {LM} becomes the GS 
of H + Zlm- The commutator [H, Zlm] vanishes indeed, 
and given A, t, v and u, there are always positive, large 
enough values for B and C that reshuffle the spectrum 
such that the lowest {LM} eigenstate ^ lm becomes the 
GS of H + Zlm under this Zeeman-like effect. We stress 
again that Zlm changes nothing in the eigenfunctions, 
eigendensities, etc., of all our Hamiltonians if w is ro- 
tationally invariant. Furthermore, angular momentum 
numbers remain approximately valid for eigenstates of 
KiMA if w is weak, and the same numbers might still 
make sense as labels by continuity when stronger defor- 
mations occur. Then the usual ad absurdum argument 
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generates a map w\ <-> t\, where T\(r) is the form fac- 
tor of the A-multipole component of the GS density for 
KlmAj leading to an exact DF, for every {LM} lowest 
state and every even A between and 2L. A generaliza- 
tion to operators Kla/a, involving (H — E n ) 2 , is trivial. 

This note offers theorems for the existence of exact 
DFs for every excited bound state, and even narrow res- 
onances, and every set of good QNs used in nuclear, 
atomic and molecular physics. Furthermore, the densi- 
ties used as arguments of our DFs do not need to be fully 
three-dimensional ones; they can be radial form factors 
of multipole components of the states under study. Our 
existence theorems, though, suffer from the usual plague 
of the field: constructive algorithms are missing and em- 
pirical approaches will have to be designed. What is the 
corresponding (KS) theory 3]? In its usual form, the 
task of calculating the kinetic energy part of the DF is 
actually left to the solution of a Schrodinger equation, 
and this can be trivially generalized to any one-body 
part. Published studies of the KS formalism are actu- 
ally dedicated to calculations of the functional derivative, 
^V xc /Sp(r), of the exchange and correlation part of the 
DF, coming from the two-body part V of the DF. Our 
present use of modified Hamiltonians, or even squares of 
H, introduces two-body operators, but also three- and 
four-body operators. For the versions where no squares 
of H occur, see the yrast suited operator K and Section 
3, the nature of the three- and four-body terms, typically 
coming from (L.L) 2 , is not forbidding, because of obvi- 
ous factorization properties. Hence a KS theory might 
be realizable for such simplified versions. With squared 
Hamiltonians, however, a KS theory seems out of reach 
at present. A systematic analysis of solvable models on 
a basis of "modes" |l9j, however, may help to extrapo- 
late such models into practical rules. For the discussion 
of differentiability, representability and fine topological 
properties of the w- and p-spaces, we refer again to 
Up to our understanding of the topology of the varia- 
tional spaces, flat or curved [2!j, of general use in nuclear, 
atomic and molecular theory, the validity domain of our 
existence theorems is quite large. We have not used the 
time dependent formalism, although much progress has 
been made in deriving excitation energies from it |30| . 
A generalization of our arguments to finite temperatures 
seems plausible, however, and insofar as inverse temper- 
ature may be viewed as an imaginary time, a generaliza- 
tion to a time dependent theory is not excluded. 
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